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ABSTRACT 

We explore the variability of quasars in the Mgii and H/3 broad emission lines and UV/optical 
continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). 

This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic 
epochs from SDSS-RM over 6 months, containing 357 quasars with Mgii and 41 quasars with H/3. On 
longer timescales, the study is also supplemented with two-epoch data from SDSS-I/IL The SDSS- 
I/H data include an additional 2854 quasars with Mgii and 572 quasars with H/3. The Mgii emission 
line is signihcantly variable (A/// ^ 10% on ~ 100-day timescales), a necessary prerequisite for its 
use for reverberation mapping studies. The data also confirm that continuum variability increases 
with timescale and decreases with luminosity, and the continuum light curves are consistent with a 
damped random-walk model on rest-frame timescales of > 5 days. We compare the emission-line 
and continuum variability to investigate the structure of the broad-line region. Broad-line variability 
shows a shallower increase with timescale compared to the continuum emission, demonstrating that 
the broad-line transfer function is not a (S-function. H/3 is more variable than Mgii (roughly by a factor 
of ^ 1.5), suggesting different excitation mechanisms, optical depths and/or geometrical configuration 
for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM 
project have important consequences for future studies of reverberation mapping and black hole mass 
estimation of 1 < z < 2 quasars. 

Subject headings: black hole physics-galaxies: active-quasars: emission lines-quasars: general-surveys 


1. INTRODUCTION 

Aperiodic luminosity variations across the electromag- 
netic spectrum are an ubiquitous feat ure of auasarJ’l 
(for a review, see lUlrich et al.l 1199711 . The optical 
and ultraviolet (UV) continuum emission of a typical 
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We use the term “quasar” to generically refer to active galactic 
nuclei with optical broad emission lines, regardless of luminosity. 


(nonblazar) quasar vary by a few tenths of a mag¬ 
nitude on timescales from weeks to years. Theoreti¬ 
cally, the observed quasar continuum variability may 
be driven by several ki nds of complex instabilities in 
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Quasar broad emission lines arise 
doppler-broadened line emission from gas deep within 
the gravitational potential well of the supermassive 
black hole, i.e., the broad line region (BLR), that is 
photoionized by the extreme UV (EUV) accretion disk 
continuum radiation. As a result, they vary in response 
to the continuum variations after a light-travel time 
delay. The amplitude and shape of the emission-line 
response are governed by th e broad-line transfer function 
(|Blandford fc McKed Il982f) . The transfer function is 
ultimately determined by the radiative mechanism, as 
well as the structure, dynamics, and ionization state 
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of the BLR. The variability time delay enables rever- 
berat ion mapping to study the structure of the BLR 
(e.g., iGaskelll 1200911 and (with assumptions about the 
geometry a nd BLR dynam i cs) estimate the black hole 
mass le.g.. IPe terson| 19931 iLaod 1199^ iPeterson et al.l 
119981 : iKaspj et all 120001 : iPetersonI 120141 . In theory, the 
reverberation mapping technique can be performed 
using any broad emission lines that respond to the 
variability of continuum emission. In practice, however, 
reverberation mapping has been largely restricted to 
using the H/? emi ssion line in low-luminos ity systems at 
low redshift le.g.. IPeterson fc Bentd[2?K^ . 

Employing reverberation mapping with optical spec¬ 
troscopy at z > 1 is critical for our understanding of the 
mass growth of supermassive black holes, as most mass 
growth occurs at this epoch. However, optical reverber¬ 
ation mapping at z > 1 requires using rest-frame UV 
broad emission lines, such as Mgii. Mgii is ionized by 
E > lb eV photons (i.e., the ionization energy to go from 
Mg II to Mg III) , which is similar to the ionization energy 
{E > 13.6 eV) of H i (which subsequently leads to recom¬ 
bination and H/3 emission). In addition, the similarity 
between the Mgii and H/? velocity widths indicates the 
two lines are produced in a similar envi ronment and dis¬ 
tance from the central black hole (e .g.. [Shen et al.1l2008l : 
IWang et al.ll2009l : iShen fc Liull201^ . However, the valid¬ 
ity of Mg II as a black-hole mass estimator is under de¬ 
bate. The variability of Mg ii has only been observed in 
a handful of quasars, and it is not entirely clear that this 
variability can be robustly traced to coher ent reverber¬ 
ation of the ionizing EUV continuum (e.g.. iClavel et al.l 


ation ot the ionizing nUy continuum (e.g., K HMmet_alJ 

iggiHReichert et al.ll994lTrevese et al.ll2007l : lWodl2008l : 


Hrvniewicz et al.l 12014 iCackett et al.l I2015D . Moreover, 

the radiative mechanism to produce Mg ii may also differ 
from that for H/3, as the former may mostly be collision- 
ally excited while the latter is a recombination line (e.g., 
lMacAlT)in^ll972|: iNetzeil 119801) . Therefore, it is vital to 
investigate the variability of Mgii for a large sample. 
This can only be accomplished with a large multi-epoch 
broad line quasar spectroscopic survey. 

In this work, we measure the ensemble variability of 
broad emission lines and continuum emission using data 
from the Sloan Digital Sk y Survey R,everber ation Map¬ 
ping project (SDSS-RM, IShen et all l2015aD and from 


supplemental observations in SDSS-1 and SDSS-H. In 
Section [2] we describe the SDSS observations. Section [3] 
presents the basic variability of quasar light curves. In 
Section IH we present the distribution of the observed 
luminosity variability. In Section [51 we introduce the 
structure function as a tool to study quasar variability 
on different timescales. Section [6] presents the structure 
function of continuum emission, and Section [7] describes 
the structure function of Mgii and H/3. In Section |8l 
we discuss the physical implications of our results. The 
main results of this work are summarized in Section O 
We adopt a flat ACDM cosmology with Dm = 0.3 and 
/iQ = 0.7. Throughout this work, “(x)” and “x” repre¬ 
sent the arithmetic mean and the median of the variable 
X, respectively. 


2. SAMPLE DEFINITION 

In this work, we use SDSS data to study quasar 
var iability. T h e SP SS 2.5 m telescop e is d escribed 
by iGunn etldl (|2006D . lEisenstein et al.l (|2011l) give a 


technical summary of the SDSS-HI project, and the 
SDSS/BOSS spectrograph and r eduction pipel i ne ar e 
described by iBolton et al.l 1)20121) : iDawson et al.l (|2013D : 
iSmee et ^ (|2013[ ). 

We focus on the data from the SDSS-RM project, 
which is an ancillary program within SDSS-HI and 
probes the variability of quasars on rest-frame timescales 
of 1 < At < 100 days. The SDSS-I/H projects provide 
ancillary data to study the variability of quasars on rest- 
frame timescales of 100 ^ At < 1000 days. 

2.1. SDSS-RM quasars 

The SDSS-RM sample, observe d during the SDSS-HI 
BOS S survey (lEisenstein et al.l 120111 : IDawson et al.l 
IMl), consists of 849 quasars, each with 32 epochs 
of observations (with a total exp osure time o f ^ 6 0 
hours): for technical details, see IShen et HI l)2015aD . 
Three out of the 32 epochs have low S/N spectra (i.e., 
S/N < 0.7 (S/N)) and are rejected. The spectrograph 
has a wavelength range of 3650-10400 A wi th a spec¬ 
tral resolution of i? ^ 2000 (jSmee et al.ll2013l ). The flux 
calibration was performed based on 70 standard stars at 
each epoch, using an improved vers ion of the standar d 
BOSS pipeline (for more details, see lShen et al.ll^l5aD . 

For each quasar, we fit the quasar spectra and obtained 
(depending on the observed spectral coverage) L = XLx 
for rest-framJn A = 3000 A and/or A = 5100 A, as well 
as the luminosities and the full-width-half-max velocity 
(ufwhm) of Mg II and/or H/3. The details of the contin¬ 
uum and the line fitting are desc ribed in detail in earlier 
work f see IShen et al.ll2008l l201ll ). 

We then selected a parent sample of 731 quasars with 
broad Mgii or H/3 by requiring z < 2.462. This require¬ 
ment ensured that either Mg ii or H/3 (or both) and their 
respective continuum regions were present in the BOSS 
spectrum. We further required that the median xfwhm 
over 29 epochs wfwhm > 1000 km s“^. 

The observed variability is a superposition of the in¬ 
trinsic variability of quasars, the measurement errors, 
and the spectrophotometric errors (see Appendix IbI). We 
applied the following sample-selection criteria to obtain 
an unbiased measurement of the intrinsic variability of 
quasars: 

• There are “dropped” epochs where the fiber was 
not properly plugged into the spectroscopic mask 
resulting in sudde n, unusually large reduction in 
the spectral flux (jShen et al.l 12015^ . To avoid 
these “dropped” spectra, we reject epochs with 
|m — to| > 1, where m is the magnitude at a given 
epoch and rh is the median magnitude over 29 
epochs. About 1% of th e total epochs are e xcluded 
1 roughly consistent with lShen et aLll2ni5aD . Visual 
inspection shows that this criterion does not reject 
any cases of real variability. 

• There are spectra with low S/N in the r-band and 
therefore with poor flux calibration, since the spec¬ 
tra are calibrated only in r-band. To identify these 
quasars, we reject quasars with r-band S/N less 
than 50, measured by convolving the spectra with 
the SDSS r-band filter. 

Throughout this work, the wavelengths of quasar features are 
always rest-frame, unless otherwise specified. 
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• There are spectra that have strong broad- 
absorption features in their Mgii profiles. In 
this case, the measurements of Mgii flux can be 
problematic. To avoid this issue, we reject these 
quasars. 

• The observed-frame quasar spectra are dominated 
by sky lines at wavelengths larger than 9000 A. 
To avoid this issue, we reject quasars with redshift 
z > 2.0 for the 3000 A continuum and Mgii and 
z > 0.8 for the 5100 A continuum and H,8. 

• There are quasars whose observed variability is 
dominated by measurement errors rather than in¬ 
trinsic variability. In this case, the estimated in¬ 
trinsic variability can be highly biased. To avoid 
this issue, we reject quasars with median signal-to- 
noise ratio (S/N) of the continuum or line luminos¬ 
ity < 10, i.e., S/N < 10. 

• There are two sources with nearby (angular dis¬ 
tance < 5 , using the SDSS imaging) bright fore¬ 
ground object that contaminate the quasar spectra 
in some epochs and introduce artificial flux varia¬ 
tions in the fibers. To avoid this issue, we reject 
these two sources. 


The final sample (i.e., clean sample) that passed the se¬ 
lection criteria and will be used for subsequent variability 
analysis of each continuum and emission-line component 
is summarized as follows: 

• the 3000 A continuum: 577 quasars; 

• Mg II broad emission line: 357 quasars; 

• the 5100 A continuum: 97 quasars; 

• H/3 broad emission line: 41 quasars. 


For each quasar in the parent sample, we estimated 
LboI using T 3000 and/or Tsioo- The bolometric correc¬ 
tion factor is assumed to b e 5 for L 3000 and 10 for T 5100 
(e.g.. [Richards et ^l2006rF^ . If a spectrum covered both 
3000 A and 5100 A we adopted the L 3000 estimator which 
is less contaminated by host galaxy starlight and there¬ 
fore provides a less biased measure of quasar luminosity. 

We also measured Mbh using the single-epoch broad¬ 
line (Mg II and/or H/1) estimators: 


, , -14bh X 3,0, , 'Cfwhm X 

( 1 ) 

For H/3, A = 5100 A, A = 6.91, B = 0.5 
(|Vestorgaard fc PetersonI 1200611: for Mgi i, A = 3000 A, 

A = 6.74. B = 0.62 (IShen et al.l[201lh . The quanti¬ 
ties AL\ and ufwhm are measured from the averages of 


Monochromatic bolometr ic correction factors are likely to 
be lu minosity-dependent (e.g.. [Lusso et ^120121 : IKrawczvk et alj 
120131) . In this work, we use Z/Boi only to divide the sources into 
different lumi nosity bins, and so we adopt the constant bolometric 
corrections of I Richards et ahl 1120061) for simplicity. Our conclusions 
do not change if we instead adopt luminosity-dependent bolomet¬ 
ric correction factors, which merely change the bin divisions by a 
small amount. 
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Fig. 1.— The distribution of the parent sample of SDSS-RM 
sources in the I/Bo1~-^bh plane. The solid, dashed and dotted con¬ 
tours in each panel show boundaries of regions containing 68.2%, 
95.4% and 99.7% of the total number of sources in our parent SDSS- 
RM sample. Points in each panel represent the SDSS-RM quasars 
that passed the selection criteria (i.e., the clean sample). The red 
dashed line in each panel represents the Eddington luminosity as 
a function of Mbh- The error bar in each panel indicates the typi¬ 
cal uncertainty of quasar luminosity (due to the uncertainty of the 
bolometric correction) and Mbh- 

29 epochs. The measurement uncertainties of and 
i^FWHM are typically 0.01 dex and 5%, respectively. The 
uncertainty of Mbh is therefore dominated by the intrin¬ 
sic uncertainty of t he estimator , which is ^ 0.4 dex (for 
a recent review, see lShenll2013ll . If a spectrum included 
both Mg II and H/3, we adopted the average Mbh (there 
are 134 such spectra i n the parent samples), following 
IVestergaard fc Osmed (|2009ll . For these sources, the dif¬ 
ference of the two Mbh estimators has a median value of 
—0.12 dex (i.e., on average, Mgii estimators give slightly 
higher Mbh than those using H/3) and a standard devi¬ 
ation of 0.35 dex. Only for three sources in the parent 
samples (but not in the clean samples) are the differ¬ 
ences in the two Mbh estimators larger than 1 dex. Two 
are z < 0.4 sources and therefore upwHM of Mgii is not 
well constrained; the redshift of the remaining source is 
z = 0.923 and therefore ufwhm of H/3 is poorly mea¬ 
sured. 

Figure [1] presents the distributions of the parent sam¬ 
ple and the clean sample in the LboI — Mbh plane for 
each component. For the 3000 A continuum and Mgii, 
the clean samples and the parent samples share similar 
parameter space. However, the clean 5100 A continuum 
and H/3 samples cover only the low-luminosity and small 
Mbh portion of parameter space, compared to the par¬ 
ent samples. This result is due to the fact that we can 
only measure the 5100 A continuum and H/3 for z < 0.8 
sources and low-redshift sources are more likely to be less 
luminous, on average, for a flux-limited survey. 


2.2. SDSS-I/II quasars 

We us e ancillary data compiled from the SDSS-I/H 
surveys (|York et al.l 120001) to measure the variability of 
quasars on longer timescales. This sample includes every 
quasar with multiple spectroscopic epochs (at least two) 
in SDSS-I/H. Most of these SDSS-I/II quasars were only 
observed twice. Only a small fraction of these quasars 
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are observed with SDSS-III. We do not supplement this 
set of quasars with SDSS-III observations due to the dif- 
ference in flux calibr ation from SDSS-I/II to SDSS-III 
(jMargala et alJ [2?)T^ . We fit the spectra and obtained 
the luminosities at A = 3000 A and A = 5100 A and 
also the luminosities of the emission lines, Mg ii an d H/3, 
following our earlier work (|Shen et alJ 1200811201 If) . For 
the 3000 A continuum and Mgii, a parent sample con¬ 
sisting of 4599 quasars was compiled. For the 5100 A 
continuum and H/3, a parent sample with 1347 quasars 
was constructed. 

To ensure that the intrinsic variability is accurately 
measured, we applied the fourth and fifth selection cri¬ 
teria described in Section [Q to our SDSS-I/II sources. 
It is not necessary to remove dropped spectra or con¬ 
strain the r-band S/N, since these problems are less fre¬ 
quent in SDSS-I/II spectra (since the spectroscopic flux 
limit is much shallower). Even if these problems occur 
in a epoch for a quasar, it only affects a single flux pair 
(i.e., / 2 // 1 ) in contrast to SDSS-RM, in which it affects 
28 flux pairs, and we use hundreds to thousands of flux 
pairs when calculating the variability. The final sample 
(i.e., clean sample) that passed the selection criteria and 
will be used for subsequent variability analysis of each 
continuum and emission-line component is summarized 
as follows: 

• the 3000 A continuum: 4213 quasars; 

• Mg II broad emission line: 2844 quasars; 

• the 5100 A continuum: 1064 quasars; 

• H/3 broad emission line: 572 quasars. 

For each quasar in the parent samples, we calculated 
Mb h an d LboI following the methods described in Sec- 
tion l2.ll Figure [2] shows the distribution of the SDSSI/H 
sources in the TboI ~ Mbh plane. Three contours again 
indicate the area that contains 68.27%, 95.45% and 
99.73% of the parent samples. The points display the dis¬ 
tributions of the final samples. For each component, the 
final sample covers similar parameter space as the par¬ 
ent SDSS-I/II sample. For the 3000 A continuum and 
Mg II, compared to the distribution of final SDSS-RM 
sources considered, most of the final SDSS-I/II sources 
cover a similar parameter space. For the 5100 A contin¬ 
uum and H/3, compared to the distribution of the final 
SDSS-RM sources, the Hnal SDSS-I/II sources cover the 
high-luminosity and large Mbh parameter space. Our 
analysis treats the SDSS-RM and SDSS-I/II datasets in¬ 
dependently, rather than combining the two datasets. 
This approach was adopted also because the two datasets 
have different time resolution and the variability of broad 
emission lines may be correlated with additional param¬ 
eters beyond LboI and Mbh- 

3. LIGHT CURVE STUDY 
3.1. Observed Light Curves 

We start with the basic properties of the observed light 
curves. Figure|3]plots continuum and broad emission-line 
light curves of the clean SDSS-RM samples. Through¬ 
out this work, we use magnitude (rather than flux or 
luminosity) changes to characterize variability so that it 
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Fig. 2.— The distribution of our SDSS-I/II sources in the 
LboI — -^BH plane. The solid, dashed and dotted contours in 
each panel show boundaries of regions containing 68.2%, 95.4% 
and 99.7% of the total number of sources in our initial SDSS-I/II 
sample. Points in each panel represent the distribution of our final 
SDSS-I/II sample after applying our selection criteria. The red 
dashed line in each panel represents the Eddington luminosity as a 
function of Mbh- The error bar in each panel indicates the typical 
uncertainty of quasar luminosity and Mbh- 




Observed—frame time since MJD 56660.209 [days] 

Fig. 3. — Continuum and emission-line light curves of typical 
SDSS-RM sources (red lines/points), and the 25-75 percentile 
range of the ensemble light curves (green shaded regions). The 
observed ensemble variability is a few tenths of a magnitude in 
both the continuum and broad emission lines. 


is straightforward to compare our results with previous 
photometric studies. The light curves do indicate that 
there is variability in both continuum emission and broad 
emission lines. However, the observed variability is a su¬ 
perposition of intrinsic variability, measurement errors, 
and spectrophotometric errors. The spectrophotomet- 
ric errors of SDSS-RM for a single epoch are at least 
^ 0.04 mag, and larger at long and short wave lengt hs, 
as computed from standard stars: see Appendix lB.il 

3.2. Basic Light-Curve Variability 
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Following iSesar et all (|2007f) , we define the basic in¬ 
trinsic variability of a light curve with n epochs as, 


O'lc — 


\ 


— ^(mi - (m))2 - (cr2)) 


( 2 ) 


i=l 


where rrii is the observed magnitude at each epoch, its 
corresponding uncertainty is Ce (ce is a summation of 
measurement errors and spectrophotometric errors in 
quadrature), and (m) is the mean magnitude across all 
epochs, aic is set to zero if the operand within the 
square root is less than zero (the fraction of such sources 
is less than 10%). This is motivated by the fact that 
the variance of the quasar variability is > 0. The me¬ 
dian ratios between the observed variance and the vari¬ 
ance due to error (measurement and spectrophotomet¬ 
ric), “ ("i))^/ (ffe), is 5.6 for the 3000 A 

continuum, 2.5 for Mgii, and 2.9 for H/3. Note that we 

required S/N > 10 to ensure that, in most cases, the 
measurement errors do not dominate the observed vari¬ 
ances. 


3.2.1. The 3000 A Continuum 

We first present the variability of the 3000 A con¬ 
tinuum as a function of quasar luminosity (Figure S]). 
We tested the correlation between the light-curve vari¬ 
ability and TboI using the Spearman rank correlation 
tesiFq. The null hypothesis of this test is that there 
is no correlation between the input datasets. The cor¬ 
relation coefficient of this test is p = —0.44. The p 
value (i.e., the probability of being incorrect in reject¬ 
ing the null hypothesis) is only 5 x 10“^®. Therefore, we 
conclude that there is a significant anti-correlation be¬ 
tween the light-curve variability and quasar luminosity 
at significance level of a = 0.00 (i.e., the probability 
threshold below which the null hypothesis will be re¬ 
jected). As we show in Section [6l this anti-correlation 
is not simply due to more luminous quasars having 
smaller Ce or typically being at higher redshift (sampling 
shorter rest-frame timescales). This behavior is consis¬ 
tent with previ ous work based on broad-band photomet¬ 
ric data fe.g.. iVyden Berk et al.l 12004 Ide Vries et al.l 
l2005t iBauer et al.ll2009HMacLeod et al.ll2012tl . 

Previous work suggests that the Eddington ratio, 
which depends on both quasar luminosity and Mbh, is 
the main driver of quasar continuum variability (e.g. , 
Wilhite et al.l I2008t IBauer et al.l I2009t lAi et al.l 120101 : 
MacLeod et al. ioig). We verify this scenario by testing 
the possible correlation between the light-curve variabil¬ 
ity of the 3000 A continuum and the Eddington ratio. 
We stress that the uncertainty of the Eddington ratio is 
rather large since both quasar luminosity and Mbh are 
uncertain by a factor of ^ 3. The Spearman correla¬ 
tion coefficient of this test is p = — 0.2 and the p value is 
5 X 10“^. This anti-correlation is consistent with previous 
work. 

We also applied the Kendall rank correlation to test our data. 
The results of this test are consistent with those of the Spearman 
rank test, unless otherwise specified. 

20 Throughout this work, we adopt a = 0.01 when we perform 
statistical hypothesis tests. 
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Fig. 4.— The light-curve variability of the 3000 A continuum. 
Each point represents a quasar. The solid, dashed, and dotted lines 
represent Icr, 2a and 3cr bounds of the quasar number distribution. 
The color indicates the probability density of each point (calculated 
via kernel density estimation). The light curve varia bility decreases 
with quasar luminosity. In this figure and Figures 15191 the error 
bars indicate the median uncertainty of the light-curve variability. 
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Fig. 5.— The light-curve variability of Mgii. Each point repre¬ 
sents a quasar. The solid, dashed, and dotted lines represent Icr, 
2cr and 30- bounds of the quasar number distribution. The color 
indicates the probability density of each point. The light curve 
variability of Mgii also decreases with quasar luminosity. 


3.2.2. Mgll 

We explore the variability of the light curve of Mg ii 
as a function of quasar luminosity (Figure [5]). With the 
same Spearman correlation test, for Mg ii, p = —0.24 and 
the p value is 3 x 10“®. Therefore, again there is an anti¬ 
correlation between the light-curve variability of Mgii 
and quasar luminosity. The same Spearman correlation 
test also suggests that the light-curve variability of the 
Mg II and the Eddington ratio are anti-correlated (with 
p = —0.24, and the p value is 7 x 10“®). As we show 
in Section El these anti-correlations are not simply due 
to more luminous quasars preferentially having smaller 
(Te or lying at higher redshift (with shorter rest-frame 
timescales). 

We compare the light-curve variability of Mgii with 
that of the 3000 A continuum in the left panel of FigureEl 
The Spearman correlation test again reveals p = 0.62 and 
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the p value is 6 x 10 i.e., our data favor a significant 

positive correlation between the light curve variability of 
Mg II and that of the 3000 A continuum. This correlation 
indicates that the variability in the light curve of Mgii 
and that of the 3000 A continuum are connected. The 
simplest explanation is that Mgii varies in response to 
the 3 000 A continuum (see the cross correlation analy¬ 
sis of IShen et al.ll20153l . Furthermore, the scatter in the 
correlation (for instance, a few sources show significant 
variability in the 3000 A continuum but the variability 
of Mg II is consistent with 0) suggests that the response 
process may not be uniform in all quasars, e.g., the re¬ 
sponse of the emission line to the continuum and/or the 
structure of the BLR may differ in different quasars. 

In the right panel of Figured we plot the difference be¬ 
tween the light curve variability of the 3000 A continuum 
and that of Mg ii versus quasar luminosity. Most quasars 
vary slightly more in the 3000 A continuum than in the 
Mg II line; the median crzc.sooo ~ czc.Mgii is 0.003 mag. A 
Spearman rank correlation test between the difference in 
the 3000 A continuum and Mg ii variability and TboI indi¬ 
cates that the difference is anti-correlated with quasar lu¬ 
minosity, with p = —0.2 and a p value of 10“^; as quasar 
luminosity increases, the variability of the 3000 A con¬ 
tinuum decreases more rapidly than the Mg ii variability. 

3.2.3. W 

We then present the light-curve variability of H/3. Since 
our final sample of H/3 consists of only 41 quasars, we 
interpret our data with caution. 

Figure [7] shows the light-curve variability of H/3 as a 
function of quasar luminosity. The Spearman correlation 
test suggests that there is no significant correlation be¬ 
tween the variability and quasar luminosity (p = —0.09 
with the p value of 0.58). This result may be due to the 
small size of our sample. We test such small-sample ef¬ 
fects by randomly selecting 41 quasars from the Mgii 
clean sample and testing the correlation between the 
Mg II light-curve variability and quasar luminosity. We 
repeated this simulation 10“* times and found that ^ 15% 
of the time, the anti-correlation between Mg ii variability 
and quasar luminosity of the limited sample is as weak as 
or weaker than the Spearman correlation test for H/3 vari¬ 
ability with quasar luminosity. From this simulation, we 
conclude that the lack of the correlation between the vari¬ 
ability of H/3 and quasar luminosity is plausibly caused 
by the small sample size. We then tested the correla¬ 
tion between the light-curve variability of H/3 and the 
Eddington ratio (p = —0.33, and the p value is 0.035) 
and found the anti-correlation is not statistically signifi¬ 
cant. We performed the same simulation and found that 
the lack of correlation is again plausibly due to the small 
sample size. 

We then compare the light-curve variability of H/3 with 
that of the 3000 A continuum. We prefer the 3000 A 
continuum instead of the 5100 A continuum because 
the galaxy contamination to the 5100 A continuum is 
more significant (and galaxy contamination biases the 
observed variability, see Section 16.211 . We can only com¬ 
pare quasars for which H,5 and the 3000 A continuum 
were both observed in the BOSS spectrum (35 quasars). 
The left panel of Figure |8] presents the light curve vari¬ 


ability of H,5 as a function of that of the 3000 A contin¬ 
uum. The Spearman correlation test indicates that there 
is no significant correlation between the variability of H,5 
and that of the 3000 A continuum (p = 0.28, and the p 
value is 0.1). Once again, we test if the lack of correlation 
is caused by small sample size by testing for correlations 
in random subsets of 35 quasars from the Mgii clean 
sample. In 10^ simulations, the Mg ii and continuum 
variability are uncorrelated ^ 0.5% of the time. From 
this we conclude that the lack of correlation between H/3 
and continuum variability is not solely due to selection 
effects. Given the fact that previous reverberation map¬ 
ping work reveals that H/3 does respond to continuum 
variability, our results indicate that H/3 and Mgii vari¬ 
ability relate to continuum variability in different ways. 
For example, the intrinsic correlation between the H,5 
variability and continuum variability might be slightly 
weaker. It is also possible that the transfer function of 
H,5 increases significantly with quasar luminosity, thus 
making the H,5 variability behave opposite to continuum 
variability as a function of quasar luminosity. The differ¬ 
ent relationships of H/3 and Mg ii variability with quasar 
luminosity in turn probe the differences between the BLR 
gas responsible for each line (see Section 1831 ) . 

In the right panel of Figure [3 we plot the difference 
between the variability of the 3000 A continuum and that 
of H/3 as a function of quasar luminosity. We again tested 
the correlation between the difference of the variability 
and quasar luminosity using the Spearman correlation 
test. We found that there is an anti-correlation between 
the differences of the variability and quasar luminosity 
(p = —0.56, and the p value is 4 x 10“"^). This correlation 
makes sense given the previously-found anti-correlation 
of the 3000 A continuum variability with quasar lumi¬ 
nosity, even with the lack of correlation of H/3 variability 
with quasar luminosity. 

We performed similar comparisons between H/3 and 
Mg II. Our results are presented in Figure HI Again, 
we can only consider quasars for which Mgii and H,5 
were both observed in the BOSS spectrum (26 quasars). 
The Spearman correlation test suggests that there is no 
strong correlation between the variability of H/3 and that 
of Mg II (p = 0.36, and the p value is 0.07), and on aver¬ 
age, the variability observed in Mg ii is less than observed 
for H/3. Note that there are seven out of 26 sources that 
show significant variability in the light curve of H/3 but 
not in that of Mgii. This is probably due to the fact 
that, as Mg II is on average less variable than H/3, the 
observed light-curve variability of Mg ii is more likely to 
be dominated by measurement and spectrophotometric 
errors. The lack of correlation between Y{(3 and Mgii 
variaiblity is likely due to the small sample size and/or 
the differences between the BLR gas that produce H/3 
and Mg II (see Section 1531 ) . 

There might be a weak correlation between the differ¬ 
ence of the variability of H/3 and that of Mg ii and quasar 
luminosity as revealed by the Spearman test (p = 0.51, 
and the p value is 0.008) although the Kendall test sug¬ 
gests that we cannot rule out the no-correlation hypothe¬ 
sis (the p value is 0.02). This correlation, if indeed exist, 
can also be explained by the anti-correlation between the 
variability of Mgii and quasar luminosity (which holds 
even we only consider these 26 quasars), even with the 
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Fig. 6.— Left: The light-curve variability of Mgii against that of the 3000 A continuum. The light-curve variability of the two components 
are correlated. Right: The difference between the light-curve variability of the 3000 A continuum and that of Mgii. The difference is 
anti-correlated with quasar luminosity. 


Fig. 7.— The light-curve variability of H/3. Each point represents 
a quasar. The solid, dashed, and dotted lines represent Icr, 2u and 
3(7 bounds of the quasar number distribution. The color indicates 
the probability density of each point. The light-curve variability 
of H/3 does not demonstrably depend on quasar luminosity which 
may due to the fact that there are only 41 quasars in the sample. 

lack of correlation of variability with quasar luminos- 
ity. 


measurement errors, and spectrophotometric errors (see 
Appendix IB]) . The distribution of Am is not Gaussian 
for any component. Actually, the distributions are bet¬ 
ter described by the Laplace (i.e., double-exponential) 
distribution. This can be explained either by the fact 
that the spectrophotometric error s are not Gaussian (see 
Appendix IbI) or, as illustrated bv IMacLeod et all (|2012li 
(their Section 3.2.2), as natural results of a superposition 
of many Gaussian distributions. 

The distributions of Am for SDSS-I/II sources are 
shown in Figure [TTJ The SDSS-I/II distributions are 
less well-sampled in Am because each quasar only has 
two epochs, compared to the 29 epochs of the SDSS-RM 
sources. In addition, the increase in variability with At 
is clearer in Figure [TT] because of the larger range of At 
in the SDSS-I/II data (which span tens to thousands of 
days). Otherwise, the SDSS-I/II variability distribution 
is similar to our SDSS-RM sources, with similar expo¬ 
nential tails. 

5. THE STRUCTURE FUNCTION: VARIABILITY AS A 
FUNCTION OF TIME SEPARATION 

In this section, we introduce the structure function, a 
tool to study quasar variability as a function of rest-frame 
time separation (At) between observations. 
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4. OBSERVED DISTRIBUTION OF AM 

Previous studies have revealed that, in general, intrin¬ 
sic quasar variability is not consistent with white noise 
but is actually a red-noise process where quasars are more 
variable on longer timescales le.g.. iMacLeod et'mi201(1I 
hereafter MIO). The basic light-curve variability we pre¬ 
sented in Section [3] is averaged over different timescales. 

Before we calculate the quasar variability as a function 
of rest-frame At (i.e., the structure function), we first 
present the distribution of luminosity variability for two 
different At bins. We calculated Am = —2.5 log(/ 2 //i) 
for all luminosity pairs separated by At. The histograms 
of Am for the continuum and broad emission lines are 
presented in Figure [TOl Note that the observed his¬ 
tograms are a superposition of the intrinsic variability. 


5.1. Definitions 

The structure function, SF(At), measures the statis¬ 
tical dispersion of observations (e.g., luminosity or flux) 
separated by time intervals. At. Many non-parametric 
statistics can be used as estimators of the statistical dis¬ 
persion of the distribution of Am, e.g., the interquartile 
range (IQR), the median absolute deviation (MAD), the 
average absolute deviation (AAD), and the standard de¬ 
viation. The IQR and MAD statistics tend to be more ro¬ 
bust against outliers or tails in the distribution than the 
AAD or standard deviation statistics: see Appendix 
for more details. In this work, we adopt the IQR es¬ 
timator as it has been adopted in previous work. In 
order to account for the presence of measurement and 
spectrophotometric errors in our data, we estimate the 
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Fig. 8.— Left: The light-curve variability of H/3 against that of the 3000 A continuum. There is no significant correlation between the 
variability of the two components. Right: The difference between the light-curve variability of the 3000 A continuum and that of H/3. The 
difference is anti-correlated with quasar luminosity. 
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Fig. 9.— Left: The light-curve variability of H/3 against that of Mgii. There is no significant correlation between the variability of the 
two components. Right: The difference between the light-curve variability of H/3 and that of Mgii. The difference is weakly correlated 
with quasar luminosity. 


intrins ic statistical dispersion following iMacLeod et al.l 

SFiQR(At) = (0.74IQR(Am))2 - (3) 

where IQR(Am) is the 25% — 75% interquartile range 
of Am and de is the total uncertainty of Am (i.e., the 
summation of measurement error and spectrophotomet- 
ric error in quadrature). The constant 0.74 normalizes 
the IQR to be equivalent to the standard deviation of a 
Gaussian distribution; the IQR is smaller than the stan¬ 
dard deviation for a Laplace distribution. Appendix 
also includes discussion of alternative structure function 
estimators: AAD, MAD, and a maximum-likelihood es¬ 
timator of standard deviation. 

The structure function can also be used to charac¬ 
terize the statistical dispersion of Am for a sample of 
many quasars with the same (or close) At. That is, 


the ensemble structure functioi0 is an average of the 
structure functions o f many sources (as demonstrated by 
IMacLeod et al.l 1^0081 individual and ensemble structure 
functions result in the same statistical properties). All 
these structure function estimators involve square roots, 
and therefore must not be negative. If measurement and 
spectrophotoinetric errors are significantly larger than 
the intrinsic variability of quasars and/or the measure¬ 
ment errors are underestimated by a factor of a few, the 
ensemble structure function would be strongly biased. 
Our sample’s S/N > 10 requirement is designed to min¬ 
imize this bias. 

Limited time sampling can cause spurious features 
in the observed structure function for individual objects 
(lEmmanoulopoulos et aTI 1201011 , but by averaging many quasar 
light-curves the ensemble structure function produces a better- 
sampled rest-frame time coverage. In Sections |6] & [3 we also 
perform simulations to account for possible sampling issues. 
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Fig. 10.— The histograms of observed Am of the SDSS-RM 
sources. In this figure, Am = —2.51og(/2//i), where /i and /2 are 
two measurements of flux separated by rest-frame At. The blue 
and red lines represent the variability of At < 25 days (the median 
At ~ 12 days for L 3000 and Mgii and 14 days for L 5100 and 
H/3) and At > 25 days (the median At ~ 40 days for L 3000 and 
MgII and ~ 50 days for I /5100 and H/3), respectively, ag and ai 
in each panel represent the standard deviation of the short and 
long timescale distributions, respectively. There are exponential 
tails in the distribution of Am. This is likely due to the superpo¬ 
sition of the intrinsic quasar variability (which is a superposition 
of many Gaussian distributions), measurement errors, and spec- 
trophotometric errors. 
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Fig. 11.— The histograms of Am of the SDSS-I/II sources. The 
blue and red lines represent the variability of At < 200 days 
(the median At ~ 20 days) and At > 200 days (the median 
At ~ 300 days), respectively. There are exponential tails in the 
distribution of Am. Also, the variability increases with At. 

5.2. The Damped Random-Walk Model 

As pointed out by previ o us work (e.g., H^llv et al.l 
120091 iMacLeod et al.l l2012t IZu et all 120131) . the light 
curves of quasar continuum emission can be well de¬ 
scribed by the DRW model. The DRW model describes 
a random process that is characterized by the following 
covariance matrix: 

C{x{U), x{tj)) = SF^ exp(-AV'r)/2 (4) 

where At = \ti — tj\. SFqo and r are the asymptotic driv¬ 
ing amplitude and the damping timescale. The structure 


function is given by 

SF{At\T, SFoo) = SFoo\/l - exp(-At/r) (5) 

At At <C T (with r typically on the order of hundreds of 
days), Eq. [5]can be simplified as, 

SF(At) = CT^At (6) 

where a = SFoo/-\/tj and the structure function at At = 
100 days, SFioo = lOd. The ensemble structure function 
for given luminosity, A, bin is given by 

pAtmax r r 

SFEN(At) = / / dAtodrdSFoa SF{Ato\T,SFoo) 

^ At min ^ ^ 

For a quasar, r and SFoo (or d) are determined by Mbh, 
A, and wavelength (e.g., MIO). For a set of quasars 
binned in a narrow range of Mbh, A, and wavelength, 
r and SFoo (or ct, SFioo) are constant and the integrals 
above reduce to a simple summation. 

6. THE STRUCTURE FUNCTIONS OF CONTINUUM 
EMISSION 

The main purpose of this work is to investigate the 
variability of broad emission lines and their connection to 
the variability of continuum emission. Therefore, we first 
present the structure function of the continuum emission. 


6.1. The 3000 A Continuum 

In this section, we present the structure function 
of the 3000 A continuum for the SDSS-RM sources. 
We first divided our sample into 4 sub-samples by 
Aboi: AboI < 10^®-^ erg s“A 10^®-^ erg s“^ < AboI < 
1045.7 ergs-1; lO^^'^ erg s-^ < AboI < 10^®-i erg s"!; 
AboI > 10"^®'^ erg s“^. The luminosity bins are con¬ 
structed to have (roughly) equal numbers of quasars in 
each. We calculated the ensemble structure function in 
each bin using the method described in Section [5T] (sub¬ 
tracting both measurement and spectrophotometric er¬ 
rors) . The binned IQR structure functions for the 3000 A 
continuum of the SDSS-RM quasars are shown by the 
solid lines in the left and center panels of Figure fT^ 
Using a Stripe 82 quasar sample, MIO found that the 
DRW parameters are correlated with quasar properties 
in the following way, 

\ogq = Qi-fQ 2 log(A^)+Q 3 (Mi-h 23 )+g 4 log( 

( 8 ) ' 

where for q = SF^c, Qi = —0.51, Q 2 = —0.479, Q 3 = 
0.131, and Q4 = 0.18; for q = t, Qi = 2.4, Q2 = 0.17, 
Qs = 0.03, and Q 4 = 0.21 (here r is in units of days). 
The uncertainty of each coefficient can be found in MIO. 
T he intrinsic scatt e rs of t he MIO relations are presented 
in IMacLeod et al.l (l2012f) . Mi is the K-corrected rest- 
frame i-band absolute magnitude. 

To compare our results with the MIO relations (i.e., 
Eq. [5]), we performed a simulation based on the MIO 
relations (hereafter “DRW simul ation”). The simulatio n 
procedures are similar to those of iMacLeod et al.l (l2012tl : 

1. We calculate a model structure function for each 
quasar using its Mbh and AboI (translated into 












































10 


rest-frame Mi assuming the i-band bolometric cor¬ 
rection factor is 12), following Eq. [8] During 
the calculation, Mbh, ^BoI, and the coefficients 
in Eq. [8] are perturbed by Gaussian noise accord¬ 
ing to their uncertainties. The predicted u (i.e., 
SFoo/'Jt) and r (perturbed by Gaussian noise ac¬ 
cording to the intrinsic scatter in Eq. H]) are trans¬ 
lated into a structure function value at each At 
using Eq. [S] 

2. In each bin of quasar luminosity at each At, we 
randomly generated a distribution of Am using 
the model structure function for each quasar in 
that bin. We then add the measurement and 
spectrophotometric errors (following Lapla ce, not 
Gaussian, distributions; see Appendix IB. II) to cre¬ 
ate the “observed” variability of the model. 

3. We calculated the structure function (using the 
IQR estimator) from the model distribution of Am, 
subtracting the median errors in quadrature in each 
bin (see Section [5T]) . 

By doing this, the simulation includes the same numeri¬ 
cal and/or time-sampling issues as the structure function 
calculated from the real data. We repeated this simula¬ 
tion 256 times to obtain a simulated structure function 
for each luminosity bin. 

Our results are shown in the left panel of Figure[T2](for 
other structure-function estimators, see Appendix |A|. It 
is evident that the variability amplitude decreases with 
quasar luminosity and increases with At. Also, the shape 
of the observed structure function in each luminosity bin 
can be reproduced well by a simulated DRW model “ob¬ 
served” like the data (i.e., with the same time sampling 
and measurement/spectrophotometric errors). To verify 
this, we performed a statistical hypothesis test. 

The null hypothesis of our test is the following; for 
each luminosity bin, the observed structure function and 
the simulated structure function share the same shape 
but with different normalization factors (in logSF(At) — 
log At space). We calculated the difference between the 
observed logSF(At) and the simulated logSF(At) and 
their uncertainties for each luminosity bin. Note that we 
only con sidered data points with At > 5 days (see Sec¬ 
tion [8T]). If our null hypothesis is true, then the differ¬ 
ence we obtained should be, statistically, a constant (i.e., 
independent of At). We adopted the chi-squared test to 
assess this hypothesis. For three of the four luminosity 
bins considered, our data fail to reject the null hypothe¬ 
sis that the two structure functions share the same shape 
(i.e., the p value of the chi-squared test is 0.01). The 
simulated and observed structure functions have different 
shapes only in the 10"^^'^ erg s“^ < LboI < 10 ^^'^ erg s“^ 
bin (although this happens only 25% of the time), but 
differ in slope by only A/3 ~ 0.1 (parameterizing each by 
SF(At) - Atf^). 

To characterize the sensitivity of our statistical test, 
we simulated two pairs of structure functions (both fol¬ 
low SF(At) ^ At^) with different slopes /3 and with the 
same S/N as our observed structure functions and our 
“DRW simulation” structure functions. We then applied 
our statistical test to each pair of simulated structure 
functions and calculated the p value of the null hypothe¬ 


sis that the shapes are the same. We repeated this simu¬ 
lation 10^ times and found that in 80% of the simulations 
(which is the most widely adopted value in statistical 
power analysis), our null hypothesis can be rejected if 
the difference in /3 is > 0.1 (> 0.2, for the highest lumi¬ 
nosity bin): the difference in slopes is < 0.1 1^3 Therefore, 
we conclude that, for each luminosity bin, the observed 
structure function and the structure function from our 
“DRW simulation” share the same shape (the difference 
in slo pe /3 < 0.1) . Tha. t is, as revealed by previous work 
(e.g.. iKellv et al.ll2009t [MacLeod et al.ll2010L l2012tl . the 
variability of the quasar continuum emission can be de¬ 
scribed by the DRW model. 

As for the variability amplitude, we find that the 
MIO relations under-predict the amplitude (by a fac¬ 
tor of ~ 1.3) except for the highest luminosity bin (i.e., 
logLBoi > 46.1). This indicates that while the MIO re¬ 
lations are accurate for high-luminosity quasars, a small 
revision is required to reproduce the variability of lower 
luminosity quasars. Note that the quasar sample used 
in MIO consists of luminous quasars {Mi < —23) while a 
significant fraction of quasars in our sample are less lumi¬ 
nous {Mi > —23). That is, our lower-luminosity quasars 
require extrapolation of the MIO relations, and this ex¬ 
trapolation seems to be inaccurate by a small factor. 

In the middle panel of Figure [121 we illustrate the 
importance of “observing” the model structure func¬ 
tion, comparing the “observed” model structure func¬ 
tion to the basic model without bias treatment. On 
long timescales (^ 20 — 100 days), the two versions of 
the structure functions are similar. On short timescales 
(< 20 days), the two versions of structure function are 
different and the differences increase with decreasing At. 
This is because, on short timescales, the uncertainties of 
our data are comparable to or even dominate over the in¬ 
trinsic variability, and the bias of the observed structure 
function is significant. In this case, if we add uncertain¬ 
ties and then subtract them via Eq|3l the model becomes 
biased toward higher amplitude variability. 

We then studied the structure function of the 3000 A 
continuum on longer {At ~ 1000 days) timescales using 
the ancillary SDSS-I/II data. Similar to our SDSS-RM 
sources, we divided the sample into four bins of LboI and 
calculated the structure function in each bin. We then 
test whether the MIO relations can effectively describe 
our data by performing the “DRW simulation” for SDSS- 
I/II sources. 

The right panel of Figure [12] shows our results. Simi¬ 
lar to the SDSS-RM sources, the structure function of 
the 3000 A continuum for SDSS-I/II quasars also in¬ 
creases with At and decreases with LboI- On timescales 
of ~ 10 — 100 days, the SDSS-RM and SDSS-I/II quasars 
have similar variability amplitudes and the MIO rela¬ 
tions again under-predict the variability amplitude by 
a factor of 1.2 for the three low-luminosity bins. The 
MIO relations successfully reproduce the shape of the 
observed structure function (the same test reveals a p 
value 0.01, and the difference in slope /3 is < 0.15) 
for each luminosity bin. In addition, using the MIO 
relation (Eq. [5|), we calculate the ensemble character- 

Throughout this work, the maximum “allowed” difference in 
/3 is estimated in this way. 
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Fig. 12.— Left panel: The structure function (solid lines) of the 3000 A continuum of the SDSS-RM sources for each lumi nosit y bin. 
The dashed lines represent our simulation results based on the MIO relations “observed” like our data (see the text in Section EUl. The 
structure functions increase with At and decrease with TboI- The shapes of the observed structure functions (on timescales At > 5 days) 
are consistent with the DRW model and the MIO relations (the difference in slope /3 < 0.1). The MIO relations slightly under-predict (by 
a factor of 1.3) the variability amplitude. Middle Panel: The dotted lines represent the structure functions from the MIO relations without 
considering the observational bias. On short timescales, the measurement and spectrophotometric uncertainties of our data dominate over 
the intrinsic quasar variability, and the structure function estimations are biased. The bias is not significant on longer timescales where 
the intrinsic quasar variability is larger. Right panel: The structure functions (solid lines) of the 3000 A continuum for SDSS-I/II sources. 
The dashed lines represent the structure functions predicted by the MIO relations. The MIO relations can roughly reproduce the observed 
structure functions. According to the MIO relations, for the four luminosity bins, the characteristic timescale r 200 days. In this figure 
and the rest of the figures, the (Icr) uncertainties of the structure functions are estimated via bootstrapping, and L represents TboI- 


istic timescale r and find that r ^ davi^ 

(where the error comes from the propagated uncertain¬ 
ties in Mbh, AboI, and the coefficients in Eq. [8]and the 
intrinsic scatter for r in Eq. [8|) . 

6 .2. The 5100 A Continuum 

We discuss the structure function of the 5100 A 
continuum for SDSS-RM quasars. Since now we are 
dealing with z < 0.8 quasars, we slightly changed 
the quasar luminosity bins: LboI < 10 "^^ ® erg s“^; 
1044.8 gj.g g-i < 2 .B 0 I < erg s-i; lO^® ^ erg s"! < 

AboI < 10"^®-® erg s“^; LboI > 10^® ® erg s“^. We calcu¬ 
lated the structure function for each bin. Our results are 
shown in the left panel of Figure [T31 As expected, the 
structure functions increase with At. In contrast to the 
structure functions of the 3000 A continuum, the struc¬ 
ture functions of the 5100 A continuum do not appear to 
be a monotonic function of quasar luminosity. 

Interpreting the structure function of the 5100 A con¬ 
tinuum is challenging because the host-galaxy contribu¬ 
tion cannot be neglected. The host-galaxy emission is 
not expected to vary on these timescales. Therefore, a 
significant host contribution would act to dilute the mea¬ 
sured Am. On the other hand, if the host-galaxy emis¬ 
sion is extended with respect to the spectroscopic fiber 
(which has a diameter of 2 ), seeing variations will cause 
apparent variability of the host-galaxy s tellar light and 
increase the observed Am (see Appendix IB.21) . 

We first restrict the sample to point sources, as clas¬ 
sified by the SDSS imaging, to eliminate luminous host 
galaxies. Restricting to point sources avoids the addi¬ 
tional variability due to variable seeing. We then re¬ 
stricted the sample to quasars with spectral decompo- 

Throughout this work, the reported uncertainties are Itr. 


sition performed bv iShen et al.l l)2015bfl and subtracting 
the estimated host contribution to the 5100 A continuum 
(quasars with > 50% galaxy light are rejected). Our re¬ 
sults are shown in the right panel of Figure [T51 

Comparing the host-corrected and uncorrected 5100 A 
continuum structure functions, it is evident that the 
galaxy emission does have significant effects: “host- 
subtracted” structure functions increase due to the re¬ 
moval of host dilution. Also, after applying the correc¬ 
tion, the structure functio ns tend to decreas e with LboI- 
Some previous works fe.g.. lKellv et al.ir2009D did not ac¬ 
count for host-galaxy emission, and therefore may have 
underestimated quasar variability at longer wavelengths. 
However, the effects of host-galaxy contamination are 
likely to be much smaller for very luminous SDSS quasars 
(e.g., MIO). 

7. STRUCTURE FUNCTIONS OF EMISSION LINES 

In Section [51 we found that the variability of the con¬ 
tinuum emission (at 3000 A) is consistent with the DRW 
model. In this section, we present the structure func¬ 
tions of broad Mgii and H/3. The study of the variabil¬ 
ity of broad emission lines provides information about 
the transfer function governing the response of the broad 
line region to incident continuum emission. The trans¬ 
fer function, in turn, reflects the structure and ionization 
state of the BLR. 

7.1. Mgii 

We first consider the variability of Mg Ii for SDSS-RM 
sources. We again divided our sources into the same four 
LboI bins as we did for the 3000 A continuum and cal¬ 
culated the structure function for each bin. We binned 
quasars by the continuum luminosity rather than the line 
luminosity so that it is easier to compare to the struc- 
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Fig. 13. — Left panel: The structure functions of the 5100 A continuum for the SDSS-RM sample. The structure functions roughly 
increase with At but are not monotonic functions of LboI- This is likely due to the contamination by galaxy emission and the effects 
of variable seeing for extended sources. Right panel: The structure functio ns of the 5100 A c ontinuum for the SDSS-RM point-source 
sample, with galaxy emission subtracted using the spectral decomposition of IShen et al.l Il2015bl ') (and removing sources with > 50% host 
contribution). After the host-galaxy effects are removed, sources in low-luminosity bins show higher variability. 


ture function of continuum emission. Our results are 
presented in the left panel of Figure [HI 

Compared to the 3000 A continuum, the structure 
functions of Mg ii also increase with At and (weakly) de¬ 
crease with LboI- We argue that this similarity is consis¬ 
tent with the idea that the variability of Mg ii is driven by 
the variability of the continuum emissio n (see the cross 
correlation analysis of iShen et al.ll201^ . 

The structure functions of Mg ii also differ from those 
of the 3000 A continuum in several aspects. First, the 
variability amplitude of the 3000 A continuum is gener¬ 
ally larger than that of Mgii, for each luminosity bin. 
In addition, the difference between the variability ampli¬ 
tudes of the 3000 A continuum and those of Mgii de¬ 
crease with LboI (similar to the results of Section [3]). 
We also found the shapes of the structure functions of 
Mg II and those of the 3000 A continuum are different by 
performing the same statistical hypothesis test we did 
in Section EH We found that, for all four luminosity 
bins, the null hypothesis is rejected (i.e., the p value of 
the chi-squared test is <C 0.01). That is, our data indi¬ 
cate different slopes for the Mgii and the 3000 A con¬ 
tinuum structure functions. Note that Mgii typically 
has larger measurement errors, and we expect it suffers 
from a stronger bias than the 3000 A continuum. To 
ensure the observed differences in shapes are not caused 
by the additional bias, we performed the “DRW simula¬ 
tion” described in Section lOI for Mgii. We found that 
this simulation cannot reproduce the observed structure 
functions of Mgii (both the normalization factors and 
shapes). Therefore, we conclude that the differences be¬ 
tween the slopes for the Mg ii and the 3000 A continuum 
structure functions are intrinsic. 


The right panel of Figure [HI shows the Mgii structure 
functions for the SDSS-I/II quasars. The variability am¬ 
plitude of Mg II is again smaller than that of the 3000 A 
continuum for each luminosity bin. We also compared 
the shapes of the structure functions of Mg ii and those 
of the 3000 A continuum by performing the same hy¬ 
pothesis test. We found that, for the two low-luminosity 
bins, the shapes of the structure functions of Mgii and 
those of the 3000 A continuum are different (i.e., the p 
value is <C 0.01). For the highest luminosity bin (i.e., the 
logLuoi >46.1 bin), the p value under the null hypothe¬ 
sis is 0.012, and we cannot reject the idea that the shapes 
are the same. If we restrict the comparison to rest-frame 
timescales At > 100 days (the timescales that are not 
covered by the SDSS-RM data), we found that we cannot 
reject the null hypothesis and the Mgii and the 3000 A 
continuum structure functions have consistent shapes in 
all luminosity bins (the p value under the null hypoth¬ 
esis is >> 0.01). The statistical test at At > 100 days 
is limited, however, by having only three data points in 
each luminosity bin and therefore lacks power to reject 
the null hypothesis (the difference in slope /3 is poorly 
constrained to be < 0.5). 

7.2. H/3 

In this section, we study the structure functions of 
H/3 for the SDSS-RM and SDSS-I/II data. In order to 
compare the structure functions of H/3 to that of the 
3000 A continuum, we only consider sources with “well- 
measured” (i.e., satisfying our selection criteria) H/3 and 
3000 A continuum. We choose the 3000 A continuum in¬ 
stead of the 5100 A continuum because the latter is sig- 
nihcantly affected by galaxy emission (see Section IQI) . 
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Fig. 14. — Left panel: The structure functions of Mg II (solid lines) for SDSS-RM sources. For comparison, we also present the structure 
functions of the 3000 A continuum (only using sources with well-measured Mgll). It is evident that the structure functions of Mg II depend 
on LboI iri the same way as the 3000 A continuum. This result is consistent with the idea that the variability of Mg II is driven by continuum 
variability. However, the 3000 A continuum is more variable than Mg II. Also, the slopes of the Mg II and 3000 A continuum structure 
functions are different (the structure function of Mg II is shallower in each luminosity bin). Note that the apparent variability decrease at 
At ~ 100 days is likely due to the fact that the spectrophotometric errors dominate over the intrinsic variability. Right panel: The structure 
functions of Mg II (solid lines) for SDSS-I/II sources. For comparison, we also present the structure functions of the 3000 A continuum 
(only using sources with well-measured Mgll). Considering the fact that SDSS-I/II sources are generally more luminous than SDSS-RM 
sources, we only created three luminosity bins. 


Due to the limited total number of sources, we do not 
divide our sources by luminosity. Figure [15] shows our 
results. 

It is evident that, for both SDSS-RM and SDSS- 
I/II sources, the variability of H/3 is smaller than that 
of the 3000 A continuum on timescales of rest-frame 
At > 20 days. We compared the shape of the structure 
function of H/3 with that of the 3000 A continuum by per¬ 
forming the same hypothesis test we did for Mgii. From 
this test, the shapes of the H/3 and continuum structure 
functions are significantly different (i.e., the p value un¬ 
der the null hypothesis is <g; 0.01). However, if we again 
restrict the comparison to rest-frame At > 100 days, the 
hypothesis that the shapes of the structure functions of 
H/3 and the 3000 A continuum are the same cannot be re¬ 
jected (the p value under the null hypothesis is 0.01). 
Again, the statistical test at At > 100 days is limited 
(the “allowed” difference in slope /3 is only constrained 
to be < 0.8). 

We also compare the structure function of Mg ii with 
that of H,5. Only 26 sources with “well-measured” Mgii 
and H/3 are included. Due to the limited sample size, 
we do not bin by luminosity. Our results are presented 
in Figure [m For the SDSS-RM sample, the variability 
amplitude of H/3 is larger than that of Mgii, which was 
also indicated by the analysis described in Section 13.21 
(Figure |9|). For the SDSS-I/II sample, the difference is 
not significant. 

We again performed the same hypothesis test to test 
whether the structure function of Mg ii and that of H/3 
share the same shape. For the SDSS-RM sample, our 
test indicates that the Mgii and H/3 structure functions 
have statistically different shapes (i.e., the p value under 



Fig. 15.— The structure functions of H/9 for SDSS-RM sources 
(pentagons) and SDSS-I/II sources (triangles). As a comparison, 
we also included the structure functions of the 3000 A continuum 
(dashed lines). H/3 is less variable than the 3000 A continuum. 
The shapes of the structure functions of H/3 and the 3000 A con¬ 
tinuum are different on timescales At < 100 days. On longer 
timescales, the shapes are similar. Note that only quasars with 
“well-measured” H/3 and 3000 A continuum are included. 

the null hypothesis is < 0.01). For the SDSS-I/II sample, 
the null hypothesis cannot be rejected by our data (i.e., 
the p value under the null hypothesis is 0.01, and the 
difference in slope (3 is 0.3). 

8 . DISCUSSION; THE NATURE OF QUASAR VARIABILITY 
8 .1. Quasar variability and quasar properties 
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Fig. 16.— A comparison of the structure functions of Mgii with 
those of H/3 for both SDSS-RM (pentagons) and SDSS-I/II (trian¬ 
gles) sources. For the SDSS-RM sources, the shapes of the struc¬ 
ture functions of Mgii and those of H/3 are different. For the 
SDSS-I/II sources, our statistical hypothesis test indicates that 
the shapes of the structure functions of Mgii and those of H/3 
are not significantly different. Note that only quasars with “well- 
measured” Mg II and H/3 are included. 


As revealed by many previous inv estiKa-tions (e.g., 
IKellv et 'n]l2009t IMacLeod et al.l[^012H . quasar variabil¬ 
ity is controlled by quasar properties (see Eq. [8|). In 
Figure na we showed that the MIO relations can effec¬ 
tively reproduce our results. Therefore, our results are 
consistent with the idea that quasar properties (e.g., LboI 
and Mbh) determine the structure of the accretion disk 
which in turn controls the instabilities in the accretion 
disk. 

The thermal timescale of the a ccretion disk is (e.g., 
iKato et ahlflMM IKellv et al.l[2n0^ 


(TH» 53.1 days, (9) 


where a is the dimens ionless viscosity parameter 
(|Shakura fc Sunvae^ll973ll . and i?s = 2 GMbh/c^ is the 
Schwarzschild radius. Assuming the 3000 A continuum 
is produced at ~ 30 (based on micr olensing obser¬ 
vations, see Eq. 4 of iMorgan et al.ir2010H . the expected 
thermal timescale is around 174 days (for a = 0.05 and 
Mbh = 10 ®'® Mq which is the median value for our 
sample). The estimated r (~ lOO.Sl^g^g^ days, see Sec¬ 
tion IHT]) from the SDSS-I/II data is consistent with the 
thermal timescale of the accretion disk. Therefore, the 
DRW model can be explained by thermal fluctu ations 
in the accretion disk (see also iKellv et al.l [200911 . On 
timescales much smaller than r, the thermal fluctuations 
(which are ultimately driven by ma gnetic turbulence , 
see recent numer ical simulations, e.g.. lHirose et al.ll2009l : 
iJiang et al.ll2013l and also a detailed theoretical calcula¬ 
tion by Lin et al. 2012) in the accretion disk drive the 
random-walk-like fluctuations in L. On longer timescales 
(^ r), the disk can adjust itself to the thermal fluctua¬ 
tions and therefore the fluctuations in L are purely white 
noise and/or are related to fluctuations in M (which cor- 
respond to the much longer visco us timescale, see e.g., 
lLvubarsidilll997t IKellv et ^1201111 . 

It has been found that, on short timescales (~ days). 


the variability of quasar continuum e mission is smaller 
than the DRW mod e l predicts (e.g.. iMushotzkv et al.l 
[20ll iZu et all IKasliwal et ahl 1201,511 . Although 

SDSS-RM includes some variability data on these short 
timescales, the median spectroscopic sampling interval 
over the campaign was ~ 4 days in observed frame (i.e., 
>1.3 days in rest frame), and the large measurement 
and spectrophotometric errors make it difficult to place 
meaningful constraints on variability over At < 5 days. 
In particular, on very short timescales (where the in¬ 
trinsic variability is small), the observed variability is 
dominated by the measurement and spectrophotometric 
errors rather than the intrinsic variability, and the struc¬ 
ture function estimation suffers from significant bias (as 
illustrated in the middle panel of Figure IT^ . 

8 .2. Implications for Reverberation-Mapping Projects 

In the SDSS-RM overview paper, iShen et ^ (l2015all 
simulated the expected results from the SDSS-RM 
project. In this section, we explore the validity of these 
simulations and the implications of our results for future 
reverberation-mapping campaigns. 

We start with the mo del used to generat e quasar con¬ 
tinuum light curves. In IShen et al.l (|2015a[l . they gener¬ 
ated continuum light curves based on the DRW model 
and the MIO relations. As shown in Section 16.11 the 
DRW model indeed describes the rest-frame 10—100 day 
variability of quasar continuum emission well. The MIO 
relations can also effectively reproduce the structure 
function, although low-luminosity quasars may actually 
be slightly more variable (by a factor of ~ 1.3) than pre¬ 
dicted by extrapolating the MIO relations. Therefore, the 
initial SDSS-RM simulations of quasar continuum light 
curves are likely to be valid. 

Justifying the transfer function is difficult since we ac¬ 
tually cannot constrain its exact shape. However, as¬ 
suming the width of the transfer function to be 10% of 
the time lag (i.e., not a S function) is consistent with 
our results. We find that the transfer function is likely 
to have a width in the range of I — 100 days (see Sec¬ 
tion [8Al). Our loose constraints on the transfer-function 
width are consistent with width estimates fr om velociW- 
resolved reverberation ma pping results le.g.. iGrier et al.l 
IM^lPe Rosa et ahllM^ . 

Last but not least, we find that, for a large population 
of quasars, Mgii varies significantly. Like H/3, Mgii is 
likely to respond to the continuum emission (as they de¬ 
pend on quasar luminosity in a similar way) and therefore 
has the potential to be used for reverberation-mapping 
campaigns. 

8.3. The Broad Emission Line Transfer Function 

The detailed ionization state, dynamical structure, and 
kinematic motions of the BLR are still poorly constrained 
from observations. In Sections [3] and [71 we compared the 
variability of two broad emission lines (Mgii and H/3) 
with that of Lsooo- We demonstrated that both broad 
lines were less variable than the continuum emission (at 
fixed time separation and luminosity). In addition, the 
structure functions of the broad emission lines are shal¬ 
lower than the 3000 A continuum structure functions, 
indicating that the broad emission lines are not driven 
by the same DRW model which describes the variability 
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of the 3000 A continuum. H/3 is typically more variable 
than Mg II. In this section, we discuss the implications 
of these results for our understanding of BLR. 

Let us first define a transfer function ^{tg) which gov¬ 
erns the response of the emission-line light curve fi {t) to 
the incident contin uum emission fr(t), after a light-travel 
time delay tg le.g.. iBlandford fc McKeelll982h : 

flit) = j dtg^itg)fcit -tg) ( 10 ) 

where /c(t) is the light curve of the continuum emission. 
The broad-line structure functions are flatter than those 
of the continuum (i.e., 3000 A, whose variability is consis¬ 
tent with the DRW model). If the variability of the EUV 
continuum is similar to the variability of the 3000 A con¬ 
tinuum, our results demonstrate that the transfer func¬ 
tions of the broad emission lines are not consistent with 
the 5-function: a broad (in time) transfer function would 
effectively flatten the input structure function. Mg ii and 
H/3 are driven by ionizing photons of E > 15 eV and 
E > 13.6 eV, respectively. This higher-energy flux likely 
has a shorter r than than the 3000 A continuum, leading 
to an apparently flatter structure function. However, ex¬ 
trapolating the MIO relations to E ~ 15 eV, the expected 
r is still ~ 100 days. 

The shape differences disappear when we only consider 
variability on timescales of rest-frame At > 100 days. 
This conclusion, if it is robust, suggests that the width 
of the transfer function for Mgii is (for our SDSS-RM 
sample) less than 100 days. The same conclusion holds 
for H/3. Therefore, on longer timescales, the variability 
of Mg II (or H/3) has similar timescale-dependence to the 
variability of the 3000 A continuum. 

Our data demonstrate that both Mgii and H/3 have 
lower variability amplitude than the 3000 A contin¬ 
uum. This result is also consistent with iMacLeod et al.l 
(see their Figure 13) who (indirectly) find that 
Mg II is less variable than the local continuum by ex¬ 
ploring the photometric va riability as a func tion of rest- 
frame wavelength (see also llvezic et al.l[2004[ l. The EUV 
continuum which actually drives both lines is proba¬ 
bly more variable than the 3000 A continuum, since 
variability increases with decreasing wavelen gth (e.g., 
IVanden Berk et all 12004 fMacLeod et al.l[2f)T^ . Indeed, 
early EUV observations reveal str ong variability (a fac¬ 
tor of 2) even wi t hin 1 day (e.g., [Marshall et al.llToOTt 
lUttlev et al.ll2000l: iHaloern et al ll2003ll . Thus the amnli- 

tude of the transfer function is likely to be significantly 
less than one, with the broad lines less variable than their 
incident continuum. This result is co nsistent with de¬ 
tailed phot oionization calculations le.g.. lKorista fc GoadI 

[IMS [200i). 

Our data suggest that the difference between the 
3000 A continuum and Mgii variability decreases with 
LboI- That is, as quasar luminosity decreases, the vari¬ 
ability amplitude of Mgii increases more slowly than 
3000 A. This may be related to a changing ionization 
structure as luminosity changes: this effectively changes 
the radius of the BLR (as inferred by reverberation¬ 
mapping studies, which measure a responsivity-weighted 
radius), and the responsivity (defined as the ratio of the 
variability of emission lines to that of continuum emis¬ 


sion! increases with increasing radius (jKorista &: GoadI 

(Mol . 

We find that H/3 is slightly more variable than Mgii, 
in agree ment with previous quasar spectral variability 
studies ([Kokubo et al.l I20I4H . The ionizing continuum 
that drives H/3 {E > 13.6 eV) is similar to that driving 
Mg II (E > 15 eV). Therefore, the differences between 
Mg II and H/3 variability are not likely to be due to the 
differences in ionizing continuum. 

The variability differences might be caused by the fact 
that H/3 is a recombination line while Mg ii is a collision- 
ally e xcited resonance line ('e.g.. lMacAlpinelll972t INetz^ 
llQSOfl . Each of these processes depends on the changes 
in incident EUV flux in a different way. This difference 
might lead to a lower variability of Mgii than H/3. 

The differences between Mgii and H/3 might also be 
caused by the different optical depth for the two emis¬ 
sion lines. Since Mgii is a resonance line, its optical 
depth is expected to be larger than that of H/3. At a suf¬ 
ficiently high density. Mg ii photons will be absorbed and 
re-emitted many times before escaping the BLR. This 
process might effectively stabilize changes in Mg ii lumi¬ 
nosity by diluting the response to continuum changes. 

The difference between Mgii and H/3 could also indi¬ 
cate the structure of the BLR that emits Mg ii is differ¬ 
ent from the BLR that produces H/3. For instance, it is 
possible the BLR material that emits Mgii extends to 
larger radius that that of H/3. Therefore, as the quasar 
continuum varies, Mgii responds on a wider range of 
timescales, and the variability of Mgii is again diluted. 

Some of these ideas are included in the “local 
optimally-emitting cloud” (LOG) model which assumes 
that the BLR consists of various density gas clouds. 
A detailed photoionization equilibrium calculation of 
the LOG model reveals that Mgii is less variable than 
H/3 over most of the typical range of BLR con ditions 
([Korista et al.lll997t iKorista &: Goadir200ClL 1200^ . 

Our data also indicate that the shapes of the structure 
functions of Mg ii and H/3 are different at At < 100 days. 
These results suggest that the transfer functions of Mg ii 
and H/3 are also different in width: that is, the transfer 
functions of Mg ii and H/3 have different radial profiles. 
At longer timescales (At > 100 days), the shapes of the 
structure functions of Mg ii and H/3 are similar (the dif¬ 
ference in slope /3 is < 0.3). This result is expected since 
the widths of the transfer functions for Mg ii and H/3 are 
both likely less than 100 days. 

Our comparison of H/3 and Mg ii variability is largely 
limited to low-luminosity quasars due to the low-redshift 
requirement for each spectrum to include both lines. 
It is possible that the relative difference in responsiv- 
ity between H/3 and Mgii has a luminosity depen dence 
(|Korista et al.l Il9^ IKorista fc GoadI [20001 1200411 . and 
comparing the different variability behavior of H/3 and 
Mg II over a larger luminosity range would place addi¬ 
tional constraints on the physical conditions and excita¬ 
tion sources of the BLR. 

9. SUMMARY 

Using SDSS-RM and SDSS I/H data, we studied the 
variability of continuum emission probed at 3000 A and 
5100 A continuum and the variability of the Mgii and 
H/3 broad emission lines. Our results can be summarized 
as follows: 











































16 


1. We determined the variability of Mgii for a large 
quasar sample. We found that, like H/3, Mgii 
varies, consistent with the scenario that Mg ii varies 
in response to the variability of the continuum 
emission (Figures 15] and [HI Sections [3] and [7]). 

2. We found that the shapes of the structure func¬ 
tions of the 3000 A continuum and those of broad 
emission lines are different, indicating the trans¬ 
fer functions governing the response of broad emis¬ 
sion lines are broad in time. We also found that 
the difference between the variability of the 3000 A 
continuum and broad emission lines decreases with 
quasar luminosity (Figures |6] and [TU Sections [3] 
and [7|) , co nsistent with photoioniz ation model pre¬ 
dictions of iKorista fc GoadI (1200411 . 

3. We confirmed that the continuum variability on 
timescales of At > 5 days is well-described by the 
DRW model (the difference in slope /3 < 0.1, see 
Figure ITT] and Section [3d]) , and that the contin¬ 
uum variability is a function of quasar properties 
(Figure [T31 Section I5TT]) . 

4. Emission from the host galaxy introduces a signifi¬ 
cant bias to the measured quasar variability at rest- 
frame 5100 A (Figure [T31 Section 13(31) . 

5. We also found that the structure functions of Mgii 
and H/3 have s tatis tica lly d ifferent shapes (Fig¬ 
ure [HI Sections 17.21 and 13(31) . Also, H/3 is slightly 
more variable than Mgii, consist ent with the pre¬ 
dictions of the LOG model (e.g.. iKorista fc Goadl 
12000L1200411 . Our results may be explained by the 
fact that the two broad emission lines have different 
radiative mechanisms, geometrical configurations, 
and/or optical depths (see Sections 131 and l8(3l) . 

Our results indicate that the predicti ons o f the SDSS- 
RM project are accurate (see Section 18.2|) . Gontinu- 
ing observations should enable accurate estimation of 
Mbh for a large set of z > 1 quasars utilizing the 
reverbera tion-mapping tech nique (using both H/3 and 
Mg II, see lShen et al.ll^l5cH . 


We thank Ghristopher Kochanek, Michael Eracleous, 
Michael Goad, Kirk Korista, and Guang Yang for ben¬ 
eficial discussions. We thank the referees for their con¬ 
structive comments. MYS acknowledges support from 
the Ghina Scholarship Gouncil (No. [2013J3009). JRT 
and YS acknowledge support from NASA through Hub¬ 
ble Fellowship grants #51330 and #51314, respectively, 
awarded by the Space Telescope Science Institute, which 
is operated by the Association of Universities for Re¬ 
search in Astronomy, Inc., for NASA under contract 
NAS 5-26555. WNB acknowledges support from NSF 
grant AST-1108604 and the V. M. Willaman Endow¬ 
ment. KDD is supported by an NSF AAPF fellowship 
awarded under NSF grant AST-1302093. 

Funding for SDSS-Hl has been provided by the Al¬ 
fred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, and the U.S. Depart¬ 
ment of Energy Office of Science. The SDSS-Hl web 
site is http://www.sdss3.org/. Funding for the SDSS 


and SDSS-H has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National 
Science Foundation, the U.S. Department of Energy, 
the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Gouncil for England. 

SDSS-Hl is managed by the Astrophysical Research 
Gonsortium for the Participating Institutions of the 
SDSS-Hl Gollaboration including the University of Ari¬ 
zona, the Brazilian Participation Group, Brookhaven 
National Laboratory, Garnegie Mellon University, Uni¬ 
versity of Florida, the French Participation Group, 
the German Participation Group, Harvard University, 
the Institute de Astrofisica de Ganarias, the Michigan 
State/Notre Dame/JINA Participation Group, Johns 
Hopkins University, Lawrence Berkeley National Labora¬ 
tory, Max Planck Institute for Astrophysics, Max Planck 
Institute for Extraterrestrial Physics, New Mexico State 
University, New York University, Ohio State University, 
Pennsylvania State University, University of Portsmouth, 
Princeton University, the Spanish Participation Group, 
University of Tokyo, University of Utah, Vanderbilt Uni¬ 
versity, University of Virginia, University of Washington, 
and Yale University. 


APPENDIX 

A. STRUCTURE FUNCTION ESTIMATORS 

The structure function measures the statistical disper¬ 
sion of two observations separated by At. Several popu¬ 
lar non-parametric statistical dispersion estimators have 
been proposed: the interquartile range (IQR), the me¬ 
dian absolute deviation (MAD), the average absolute de¬ 
viation (AAD), and the standard deviation. Their defi¬ 
nitions are: 

• the IQR es tima tor (used bv lMacLeod et al1l2012f) : 
see Section 15(11 


the MAD estimator: 


SFMAD(At) = V(1.48MAD(Am))2 - al. 


where MAD (Am) = 


Am — Am 


• the AAD estimator (used by iVanden Berk et ^ 

IHqI: 


SFAAD(At) = W-(|Am|) -(cr2). 


• the maximum-likelihood standard deviation (ML) 
estimator: the standard deviation can be measured 
using the maximum- likelihood method proposed by 
lAlmaini et al.l (j200(1f ) (see their Section 3.1). 

We subtract the spectrophotometric error from each 
structure function calculated from the stars using the 
same estimator. The constant factors in each of the 
four structure-function estimators normalizes them to 
be equivalent to the standard deviation of a Gaussian 
distribution (assuming Gaussian measurement and spec¬ 
trophotometric errors). Our data and uncertainties have 
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log(i) <45.3 



46.1>log(L) >45.7 
log(L) >46.1 



Rest—frame At [days] 


Fig. 17.— The structure function of the 3000 A continuum for 
each quasar luminosity bin. The IQR, MAD, AAD, and ML es¬ 
timators are shown. The IQR and MAD estimators are stable 
against measurements with relatively large uncertainties and also 
robust against outliers. 


non-Gaussian tails, which lead to differences in the esti¬ 
mators: for distributions with larger tails than a Gaus¬ 
sian, the IQR and MAD estimators are smaller than the 
AAD and standard deviation estimators. 

The four structure functions of the 3000 A continuum 
are plotted in Figure [iTl Note that if the number of 
flux pairs for a given At bin is less than 10, we do not 
calculate the variability in this particular bin. For each 
LboI bin with many quasars, each quasar has different in¬ 
trinsic variability and measurement/spectrophotometric 
errors. That is, the ensemble distribution of Am is a su¬ 
perposition of many Gaussian distributions and is not a 
perfect Gaussian but instead has exponential tails. More¬ 
over, the subtraction of measurement and spectrophoto- 
metric errors (i.e., subtract de or (ue)) has an effect on 
the estimation of the ensemble structure function. This 
results in slightly different results for each of the four 
structure-function estimators. 

For the AAD estimator, since we are subtracting the 
mean value of measurement and spectrophotometric er¬ 
rors, our results may be biased by quasars with mea¬ 
surement and spectrophotometric errors that are much 
larger than the intrinsic variability. The standard devia¬ 
tion estimator uses individual errors (using the maximum 
likelihood method), and so is stable against quasars with 
relatively large measurement and spectrophotometric er¬ 
rors. The IQR and MAD estimators use median values of 
the measurement and spectrophotometric errors, and so 
are also stable against quasars with large measurement 
and spectrophotometric errors. Moreover, the IQR and 
MAD estimators are also robust against outliers in the 
distribution of Am. 

The IQR, MAD and ML (or standard deviation) esti¬ 
mators give similar results: (1) the ensemble structure 
function increases with At in a similar way; (2) the en¬ 
semble structure function decreases with quasar luminos¬ 
ity in a similar way. 

The four structure functions of the 5100 A continuum 
are plotted in Figure [TH] for the full “clean” sample of 
SDSS-RM quasars. Figure [19] shows only point sources 
with galaxy emission subtracted, using the spectral de- 
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Fig. 18. — The structure function of the 5100 A continuum for 
each luminosity bin. The IQR, MAD, AAD, and ML estimators 
are shown. 
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Fig. 19. — The structure function of the 5100 A continuum for 
each luminosity bin. The IQR, MAD, AAD, and ML estimators 
are shown. Only point sources are included, and galaxy emission 
is subtracted. 


composition of iShen et al.l (l2015bll and removing quasars 
with > 50% host contribution at 5100 A. 

The four structure functions of Mg ii and H/3 are shown 
in Figures [201 and [211 respectively. In all four cases, the 
broad-line structure functions are flatter than the con¬ 
tinuum structure functions. 

As shown in this section, the IQR and MAD estimators 
are robust against o utliers. We pre f er the IQR since it 
has been adopted in iMacLeod et al.l (|2012t) . 

B. SDSS-RM SPECTROPHOTOMETRIC UNCERTAINTY 

In this section, we quantify the spectrophotometric er¬ 
rors in the SDSS-RM survey. The spectra were calibrated 
using 70 standard stars. Therefore, the variability of 
standard stars used in the SDSS-RM survey should re¬ 
flect the spectrophotometric errors. Note that one star 
was found to be intrinsically variable, and five other stars 
have an epoch with a “dropped” spectrum. We do not 
remove these stars when accounting for the variability of 
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Fig. 20. — The structure function of Mgii for each luminosity 
bin. The IQR, MAD, AAD, and ML estimators are shown. 
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Fig. 21.— The structure function of H/3 for each luminosity bin. 
The IQR, MAD, AAD, and ML estimators are shown. 


standard stars since these stars were not removed when 
calibrating the spectra. Due to their small number, these 
stars do not signihcantly affect the flux calibration. 


B.l. STRUCTURE FUNCTION OF STANDARD 
STARS 

As a first step, we present the observed distribution of 
Am for standard stars. To do so, we created 21 wave¬ 
length bins starting from 3700 A, each with a width of 
300 A. We then calculated the average flux in each bin at 
every epoch for every standard star. As a final step, we 
again calculated Am for any two observations separated 
by At. 

In Figure 1211 we plot the distribution of Am for stan¬ 
dard stars. It is evident that the variability of standard 
stars depends strongly on wavelength. The variability at 
^ 7000 A is the smallest 0.06 mag). Also, the variabil¬ 
ity increases significantly to both the short-wavelength 
0.1 mag) and the long-wavelength ends (~ 0.07 mag). 
This is due to the fact that the flux calibration is done 
using r-band only. Note that the distributions are not 
Gaussian but are instead Laplacian which indicates that 


Wavelength [A] 



Fig. 22.— The distributions of Am for standard stars. Again, 
Am = —2.51og(/2//i). There are exponential tails in the distri¬ 
bution of Am for standard stars. The variability is the smallest 
at ~ 6000 - 7000 A (since the flux calibration is done using the 
r-band), and higher at both small and large wavelengths. 

the spectrophotometric errors are not perfectly Gaussian 
(and the tails are not due to the variable stars). 

Since the quasar spectra were calibrated using stan¬ 
dard stars, we can then use the structure function of 
standard stars as an estimation of spectrophotometric 
errors. Note that the structure function quantifies the 
variability of a flux pair, and therefore the spectropho¬ 
tometric uncertainty for a single epoch is the structure 
function divided by \/2. The variability of standard stars 
depends on wavelength, and so we calculated the struc¬ 
ture function of standard stars in each wavelength bin 
separately. Figureplots our results. It is evident that 
the structure function of the standard stars is largely con¬ 
stant with At for each wavelength bin. There is a small 
increase (within 0.02 mag) at At > 20 days: this is asso¬ 
ciated with the time between different dark runs in the 
SDSS-RM observations. The structure function of stan¬ 
dard stars depends on wavelength and is the smallest in 
the ~ 6000 — 7000 A band. In the ^ 6000 — 7000 A band, 
the scale of the variability is only ~ 0.06 mag (for the 
IQR estimator). Therefore, for a single epoch, the spec¬ 
trophotometric error of our SDSS-RM data is expected 
to be not s maller than ^ 0.06 /-\/2 pc 0.04 mag (also see 
Figure 20 of iShen et al.ll20I^ . 

The spectrophotometric error of the SDSS-I/II data 
for a single epoch is assum ed to be 0.04 mag 
(jAdelman-McCarthv et al]l2008f l. 

B.2. 03 

We also checked the variability of 03A5OO7 (hereafter 
03). Physically, we expect that there is no intrinsic vari¬ 
ability of 03 on timescales of ^ 100 days. That is, the 
variability of 03 should be equivalent to that of standard 
stars. To verify this, we calculated the structure function 
of 03 while subtracting only measurement errors (i.e., 
Ue in Section l5 .1 1 only accounts for measurement errors) 
and compared the structure functions with the expected 
spectrophotometric errors of the flux pairs. The expected 
spectrophotometric errors are the average of the struc¬ 
ture function of standard stars at the same wavelength 
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Fig. 23.— The structure function of standard stars for each wave¬ 
length bin. The spectrophotometric uncertainty of a single epoch 
is the structure function divided by y/2. The spectrophotomet¬ 
ric errors depend on wavelength and have a minimum value at 
6000 — 7000 A (of 0.06/\/2 mag). Furthermore, the structure 
function only slightly depends on At (increasing by 0.02 mag at 
At > 20 days). 

as the 03 (i.e., 5007(1 + z) A, where z is the redshift 
of the quasar). We divided our sources into subsamples 
of point sources and extended sources in order to inves¬ 
tigate the additional variability due to variable seeing. 
To investigate the effects of measurement errors on the 
estimation of the structure function, we binned quasars 

by S/N. For each sub-sample, we created three bins: 
15 < ^ < 25^ < ^ < 35, ^ > 35. We required 

a minimum of S/N > 15 so that the measurement errors 
do not dominate over the spectrophotometric errors. We 
then calculated the structure function of each bin. 

Figure [M] shows the structure function for the sub¬ 
sample of point sources. The variability of [OIII] for 
point sources is nearly identical to the expected spec¬ 
trophotometric errors. For each bin, we calculated the 
difference between the structure function and the spec¬ 
trophotometric errors in quadrature. The difference is no 
more than 0.03 mag, and there is no systematic offset. 
In addition, the variability is not a monotonic function 

of S/N which suggests our S/N > 15 cut reliably pre¬ 
vents a bias from large measurement errors dominating 
the observed variability. Therefore, we conelude that our 
estimation of the speetrophotometric errors is robust. 

Figure [25] plots the structure function for the sub¬ 
sample of extended sources. It is evident that, for each 
bin, 03 shows excess variability compared to the spec¬ 
trophotometric errors. This is due to additional variabil¬ 
ity 0.07 — O.I mag in the IQR) from variable seeing 
in addition to the spectrophotometric errors. Additional 
variability from seeing changes has consequences for the 
structure-function estimates of extended sources (e.g., if 
there is a substantial host-galaxy component to the mea¬ 
sured 5100 A luminosity, see Section 15^ . 


Fig. 24.— The structure functions of 03 for point sources. Note 
that the structure functions presented here were only corrected 
for measurement errors. The dashed lines represent the expected 
spectrophotometric errors which are obtained according to the red- 
shift distribution of quasars in each bin. It is evident that, for point 
sources, 03 does not show intrinsic variability. This result indicates 
that our estimation of the spectrophotometric errors is robust. 



Fig. 25.— The structure functions of 03 for extended sources 
only. Note that the structure functions that are presented here 
were only corrected for measurement errors. The dashed lines rep¬ 
resent the expected spectrophotometric errors which are obtained 
according to the redshift distribution of quasars in each bin. Ex¬ 
tended sources show variability beyond the spectrophotometric er¬ 
rors which is likely due to seeing effects (i.e., the variable seeing in 
different epochs). 
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